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Abstract 

Describing the hyperonic and quark phases of neutron stars with an isospin- and momentum- 
dependent effective interaction for the baryon octet and the MIT bag model, respectively, and using 
the Gibbs conditions to construct the mixed phase, we study the energy release due to the hadron- 
quark phase transition. Moreover, the frequency and damping time of the first axial u>-mode of 
gravitational waves are studied for both hyperonic and hybrid stars. We find that the energy release 
is much more sensitive to the bag constant than the density dependence of the nuclear symmetry 
energy. Also, the frequency of the w-mode is found to be significantly different with or without 
the hadron-quark phase transition and depends strongly on the value of the bag constant. Effects 
of the density dependence of the nuclear symmetry energy become, however, important for large 
values of the bag constant that lead to higher hadron-quark transition densities. 
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I. INTRODUCTION 



Neutron stars (NSs) are among the most mysterious objects in the Universe. With a 
typical mass of about M = 1.4M but a radius of only about 12 km, the average density in 
NSs is several times that in atomic nuclei. Their extreme compactness makes them a natural 
laboratory to test our knowledge about general relativity and properties of dense neutron- 
rich nuclear matter Moreover, their generally large angular momentum and possible 
quadrupole deformation make them strong candidates among many possible sources to emit 
gravitational waves (GWs) 2j. It is well known that the estimated central density of some 
NSs may be higher than the predicted hadron-quark phase transition density, and thus there 
might be a quark core in these so-called hybrid stars, see, e.g., ref. jsj] for a recent review. 
Moreover, some originally hadronic NSs with central densities below but close to the hadron- 



br example, their 
5|. Therefore, the 



quark transition density might increase their central densities due to, 
spin-downs 4( or the accretion of masses from their binary companions 
hadron-quark phase transition may occur in their cores, turning them into hybrid NSs as they 
evolve. Consequently, a micro-collapse is expected to occur as a result of the softer equation 
of state (EOS) of quark matter than that of hadronic matter. Because of the difference in 
the binding energies for the pure hyperonic and hybrid configurations of NSs, some energy 
will be released from the NS after the hadron-quark phase transition. To understand the 
mechanism of the hadron-quark phase transition and the associated energy release, their 
dependence on properties of the dense hadronic and quark matter, the way the energy 
release being dissipated in NSs or carried away by gravitational waves, and the gravitational 
wave signatures of the EOS of dense matter and the expected hadron-quark phase transition 
are among the many interesting questions currently under intense investigations in neutron 
star physics and gravitational wave astronomy, see, e.g., refs. ,BJ3. 



In a recent work involving some of us 



13], several model EOSs for hybrid stars were 



obtained from an isospin- and momentum-dependent effective interaction (MDI) 13J, uA\ 
for the baryon octet, the MIT bag model for the quark matter flol and the Gibbs con- 



struction for the hadron-quark phase transition |17l . [18|. Since there have been many other 
studies using similar approaches in the literature, it is especially worth mentioning that 
while the isospin-symmetric part of the hadronic EOSs used in this work is constrained 
by the experimental data on collective flow and kaon production in relativistic heavy-ion 
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collisions 



19| . the symmetry energy term at sub-saturation densities is constrained in a 



narrow range by the experimental data on isospin diffusion in heavy-ion collisions at inter- 
mediate energies 20|. The obtained EOSs with various bag constants were then used in 
studying the main features of the hadron-quark phase transitions and the resulting mass- 



radius corre 
refs. 



13 
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ations of hybrid stars. With this paper as a continuation of our recent work in 



25] , we make new contributions to the fast growing field of nuclear astrophysics 



in two ways. First of all, critical for understanding the physics of the hadron-quark phase 
transition are the EOSs of both the hadronic and quark matter. Despite of the great efforts 
made by both the astrophysics and the nuclear physics community, our current knowledge 
about the EOS of dense matter either in the hadronic or quark phase is still very poor. More 
specifically, for the neutron-rich hadronic matter the most uncertain part of the EOS is the 
density dependence of the nuclear symmetry energy which encodes the energy associated 
with the neutron-proton asymmetry 20J. For the quark matter within the MIT bag model, 
the most uncertain part of its EOS is the bag constant which measures the inward pressure 
on the surface of the bag to balances the outward pressure originating from the Fermi motion 
and interactions of quarks confined in the bag. It is thus of interest and also necessary to 
examine the relative effects of the density dependence of the nuclear symmetry energy and 
the bag constant on the total energy release due to the hadron-quark phase transition in 
NSs. Secondly, a part of the energy release may be carried away by GWs. The maximum 
amplitude of the latter can then be obtained by assuming that all the released energy is 
available for GWs. Besides testing a fundamental prediction of general relativity, gravi- 
tational waves hold the great promises to open up a completely new non-electromagnetic 
window into the Universe. On the other hand, if detected in the future, the GWs may also 
help probe properties of dense neutron-rich nuclear matter in NSs. Thus, it is important to 
examine the imprints of the nuclear symmetry energy and the bag constant on the frequency 
and damping time of GWs. In this work, we focus on the first axial w-mode of GWs due to 
the disturbance of the space-time curvature 2a] . We find that the energy release obtained 
using the Gibbs construction for the quark-hadron phase transition is much more sensitive 
to the bag constant than the density dependence of the nuclear symmetry energy. Further- 
more, the frequency of the w-mode is found to be significantly different with or without the 
hadron-quark phase transition and depends strongly on the value of the bag constant B for 
a given density dependence of the nuclear symmetry energy. However, when a larger bag 
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constant is used such that the hadron-quark phase transition happens at a higher baryon 
density, effects of the density dependence of the nuclear symmetry energy also become very 
important. 

This article is organized as follows. In Section II, we summarize the MDI interaction, the 
MIT bag model, and the resulting model EOSs for hybrid stars obtained using the Gibbs 
construction. In Section III, the total energy release from hadron-quark phase transition 
in NSs and its dependence on the symmetry energy and the bag constant are presented. 
In Section IV, the frequency and damping time of the axial u>-mode of GWs and their 
dependence on the symmetry energy and the bag constant are studied. Finally, a summary 
is given in Section V. 



II. MODEL EQUATIONS OF STATE FOR HYBRID STARS 



Since we are mostly interested in examining the relative effects of the nuclear symmetry 
energy and the bag constant on the total energy release during the hadron-quark phase 
transition in neutron stars as well as the frequency and damping time of the u>-mode of 
emitted GWs, we first illustrate their effects on the EOS of hybrid stars by using the model 
introduced in ref. 13j. For completeness and to facilitate the discussions, we briefly sum- 



marize the EOSs for hybrid stars obtained in ref. [13[ using the MDI interaction for the 
baryon octet, the MIT bag model for the quark matter and the Gibbs construction for the 
hadron-quark phase transition. We refer the reader to ref. [l3] for details. 

Assuming that the nucleon-hyperon and hyperon-hyperon interactions have the same 
density and momentum dependence as the nucleon-nucleon interaction with the interaction 
parameters fitted to known experimental data at normal nuclear matter density, the potential 
energy density of a hypernuclear matter due to interactions between any two baryons is 
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where b (b') denotes the baryon octet included in the present study, i.e., N, A, S, and H. 
The conventions for the isospin are Tn = —1 for neutron and 1 for proton, ta = for A, 
= —1 for £~, for E° and 1 for S + , and r= = —1 for E~ and 1 for S°. The total baryon 




FIG. 1: (Color online) Density dependence of nuclear symmetry energy from the MDI interaction 
with parameter x = and —1. 



density is given by p = J2bJ2 Tb Pm an ^ fn(r,P) is the phase-space distribution function 
of particle species r&. For hyperons in symmetric nuclear matter at saturation density, 



their potentials are fixed at U^\p N = p ) = —30 MeV, U^'(p N = p ) = —18 MeV and 
Ujf (p N = p ) = 30 MeV for A, H and S, respectively. The parameter x is used to model 
the isospin dependence of the interaction between two baryons. Its value can be adjusted 
to mimic different density dependence of the nuclear symmetry energy and is taken to be 
or —1 in the present study. Shown in Fig. [T]is the corresponding density dependence of the 
nuclear symmetry energy E syra (p). It is worth noting that the latest experimental constraints 
on E sym (p) at sub-saturation densities are consistent with but span a region larger than the 
one covered by the x = and x — — 1 curves 27(. The experimental constraints on E sym (p) 
at supra-saturation densities are, on the other hand, still very uncertain. In fact, the high 
density behavior of E sym (p) is considered the most uncertain part of the EOS of dense 



neutron-rich nuclear matter 28( . One of the purposes of this work is to study effects of the 
symmetry energy. As shown in Fig. [TJ the choice of x = and x = —1 allows the variation 
of the symmetry energy in a large range at supra-saturation densities. 

To obtain the complete EOS for hybrid stars, besides the EOS for hyperons described 
above, the MIT bag model for the quark matter H] and the Gibbs construction for the 



hadron-quark phase transition 17|,ll8( have been used in ref. 



e.g., refs. 



29 



131 ] . As in previous studies, see, 



301 ]. the hybrid star is divided into three parts from the center to the surface: 



the liquid core, the inner crust, and the outer crust. It is in the liquid core where the matter 
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FIG. 2: (Color online) The EOSs of pure npefi matter (nucleonic), hyperonic matter (MDI Hyp-R 
interaction) and hybrid stars with B 1 ^ = 180 MeV and 170 MeV, and symmetry energy parameter 
x = (left window) and x = — 1 (right window). 

can be pure hadrons, quarks or a mixture of the two. In the inner crust, a parameterized EOS 
of P = a + be 4 / 3 is used as in refs. 



29 



30j. The outer crust usually consists of heavy nuclei 



and the electron gas, where the BPS EOS [3l| is used. The transition densit y a, between 



the liquid core and the inner crust has been consistently determined in refs. [29l . |30| . Taking 
the density separating the inner from the outer crust to be p out = 2.46 x 10~ 4 fm~ 3 , the 
parameters a and b can then be determined using the pressures and energy densities at 
p t and p out . Shown in Fig. [2] are the EOSs with x = and x = —1, respectively, both 
using the MIT bag constant B l / A = 180 MeV and 170 MeV. For comparisons, the pure 
npep (labeled as nucleonic) and hyperonic (labeled as MDI Hyp-R interaction) EOSs are 
also included. As it is well known, the appearance of hyperons and the hadron-quark phase 
transition softens the EOS of neutron star matter. Also, it is worth noting that the adiabatic 
coefficient 7 = d\og(P) / dlog(p) at saturation density is 2.63 and 2.57 for x = and x = —1, 
respectively. The EOSs with x = and x = — 1 are thus about the same below and around 
the saturation density. However, it is interesting to see that the starting point and the degree 
of softening due to the appearance of hyperons are sensitive to E sym (p) at high densities. 
Moreover, the E sym (p) also affects appreciably the starting point of the hadron-quark mixed 
phase especially with the larger bag constant. Nevertheless, it is obvious that the starting 
point is much more sensitive to the bag constant B for a given symmetry energy parameter 



x. These features are consistent with those first noticed by Kutschera et al. 32|. It is also 



important to emphasize that since the mixed phase is described using the Gibbs instead 



of the Maxwell construction, the energy density thus increases continuously with pressure 



across the mixed phase without a jump 



stable hybrid star must satisfy 



33[. Consequently, the Sidov criterium that a 



i.e., ^ < 



|(1 + — ) with P, e g and e# denoting 



respectively, the pressure, the energy density of quarks and hadrons at the transition point, 
is automatically satisfied. 



III. TOTAL ENERGY RELEASE DUE TO HADRON-QUARK PHASE TRANSI- 
TION IN NEUTRON STARS 



The amount of energy release during the micro-collapse of a neutron star triggered by 
the hadron-quark phase transition is given by the change in its binding energy before and 
after the phase transition. The binding energy of a neutron star can be obtained from 
first solving the Tolman-Oppenheimer-Volkoff (TOV) equation with the corresponding EOS 
(G = c = l), 

dP(r) (m(r) + 47rr 3 P(r))(e(r) + P{r)) 



dr 
dm(r) 
dr 



r(r — 2m(r)) 



(2) 



Aire(r)r 



where P(r) and e(r) are the pressure and the energy density at radius r. The binding energy 
of a neutron star is then given by [35 1 



ban 



where M g is the gravitational mass of the neutron star measured from infinity 36|], i.e., 

-R 



M n 



47rr e{r)dr. 



(3) 



(4) 



The Mb ar is the baryonic mass of the neutron star, namely, the mass when all the matter in 



351 ] . It can be calculated from Mb ar = NMa, where 



the neutron star is dispersed to infinity 
Mb = 1.66 x 10~ 24 g is the nucleon mass and N is the total number of baryons 35J, i.e., 



N 



R 



Airr' 



2m(r) 



-1/2 



p B {r)dr 



(5) 



with pb{t) being the baryon density profile of the neutron star. While the total energy 
release E g during the phase transition is the difference of Eq. ([3]) before and after the micro- 
collapse, it reduces to the difference in gravitational mass for a hadronic (h) and a hybrid 



star (q) as a result of baryon number conservation, namely, 

E g = M g>h - M m . 



(6) 
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FIG. 3: Total energy release due to the hadron-quark phase transition as a function of the baryonic 
mass of a neutron star. 

TABLE I: The baryon number (N), gravitational masses of hyperonic stars (M g (H)), hybrid stars 
(M g (HQ)), the quark core including the mixed phase (M 5 (core)), the radii of hyperonic (R(H)) 
and hybrid stars (R(HQ)) as well as the energy release (E g ) with x = and = 170 MeV. All 
masses are in unit of M & and radii in km. 



N 


Mg(H) 


Mg(HQ) 


M g (core) 


R(H) 


R(HQ) 


\og{E g ){erg) 


1.4E57 


1.0861 


1.0798 


0.886659 


11.2917 


10.1679 


52.0517 


1.5E57 


1.1570 


1.1491 


0.971092 


11.3044 


10.1256 


52.1492 


1.6E57 


1.2270 


1.2173 


1.054172 


11.3041 


10.0691 


52.2379 


1.7E57 


1.2960 


1.2843 


1.135890 


11.2909 


9.9935 


52.3199 


1.8E57 


1.3641 


1.3501 


1.217399 


11.2620 


9.8796 


52.3977 



Shown in Fig. E] is the energy release as a function of the baryonic mass Mb ar /M of 
a neutron star. To be more quantitative, shown in tables [J [Til HH an d IIVI are detailed 
comparisons of the baryon number (N), gravitational masses of hyperonic stars (M g (H)), 
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TABLE II: Same as H but with x = and B 1 ^ = 180 MeV. 



N 


M g (H) 


Mg{HQ) 


M g (core) 


R(H) 


R(HQ) 


log(E g )(erg) 


1.4E57 


1.0861 


1.0858 


0.335018 


11.2917 


11.0960 


50.7137 


1.5E57 


1.1570 


1.1564 


0.457409 


11.3044 


11.0076 


51.0370 


1.6E57 


1.2270 


1.2259 


0.582411 


11.3041 


10.8971 


51.2965 


1.7E57 


1.2960 


1.2942 


0.708860 


11.2909 


10.7625 


51.5124 


1.8E57 


1.3641 


1.3612 


0.838033 


11.2620 


10.5902 


51.7018 



TABLE III: Same as U but with x = -1 and B 1 ^ = 170 MeV. 



N 


Mg(H) 


Mg(HQ) 


M g (core) 


R(H) 


R(HQ) 


log(E g )(erg) 


1.4E57 


1.0911 


1.0813 


0.865139 


12.5570 


10.3024 


52.2408 


1.5E57 


1.1632 


1.1507 


0.953706 


12.5975 


10.2280 


52.3516 


1.6E57 


1.2346 


1.2189 


1.040034 


12.6178 


10.1446 


52.4498 


1.7E57 


1.3052 


1.2859 


1.124924 


12.6185 


10.0418 


52.5385 


1.8E57 


1.3749 


1.3516 


1.209269 


12.5998 


9.8992 


52.6205 



hybrid stars (M g (HQ)), the quark core including the mixed phase (M s (core)), the radii of 
hyperonic (R(H)) and hybrid stars (R(HQ)) as well as the energy release (E g ). It is seen 
that the energy release increases with M^/Mq and is higher with the smaller (Z? 1 / 4 = 170 
MeV) bag constant B but stiffer (x = —1) symmetry energy. Effects of varying the bag 
constant B are obviously more significant than varying the symmetry energy parameter x, 
especially on the core mass and thus the energy release. The variation of the bag constant 
B also affects appreciably the radii of hybrid stars. On the other hand, the variation of the 
symmetry energy parameter x only has appreciable effects on the radii of both hyperonic and 
hybrid stars. It's effects on the energy release is much smaller than the bag constant B. It is 
worth noting that it was first shown in ref. that the binding energies of NSs consisting 
of pure npefi matter depend strongly on the density dependence of the nuclear symmetry 
energy. It is easily understandable that the difference in binding energies before and after 
the hadron-quark phase transition becomes less sensitive to the symmetry energy. Some 
of these features can be better understood from inspecting the gravitational mass-radius 
and mass-central density correlations with and without the hadron-quark phase transition 
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TABLE IV: Same as Q] but with x = -1 and B 1 ^ = 180 MeV. 



N 


Mg(H) 


Mg(HQ) 


M 5 (core) 


R(H) 


R(HQ) 


\og(E g )(erg) 


1.4E57 


1.0911 


1.0904 


0.389260 


12.5570 


11.8321 


51.0467 


1.5E57 


1.1632 


1.1616 


0.563910 


12.5975 


11.4949 


51.4673 


1.6E57 


1.2346 


1.2314 


0.716042 


12.6178 


11.2027 


51.7548 


1.7E57 


1.3052 


1.3999 


0.855991 


12.6185 


10.9207 


51.9731 


1.8E57 


1.3749 


1.3670 


0.990235 


12.5998 


10.6140 


52.1525 



using different bag constants and symmetry energy functionals. Shown in Fig. H] are the 
gravitational mass M g as a function of radius R and central density p c for several relevant 
cases. It is seen that with B 1 ^ = 170 MeV the phase transition happens at a significantly 
lower density, resulting in a smaller maximum mass than the case with B 1 ^ = 180 MeV. 
On the other hand, for a given bag constant B, the hyperonic EOS is stiffer with a stiffer 
symmetry energy, leading to a larger radius and thus a lower central density and a smaller 
maximum mass. The mass difference before and after the hadron-quark phase transition is 
therefore also larger. 
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FIG. 4: (Color online) The mass-radius and mass-central density correlations for pure npe/j, matter 
(nucleonic), hyperonic matter (MDI Hyp-R interaction) and hybrid stars with l? 1 / 4 = 180 MeV 
and 170 MeV, and symmetry energy parameter x = (left window) and x = — 1 (right window). 

To see how the total energy release is related to emitted GWs, we add a few comments 
as follows. They are also useful for explaining why we study in detail in the next section 
the frequency and damping time of the w-mode of GWs. Generally, the magnitude of GWs 
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is denoted by the gravitational strain amplitude in the following wave form: 

h(t) = h e- {tM sm(uJot), 



(7) 



where h is the initial amplitude, Uo is the angular frequency and r gLU is the damping 
timescale. Assuming all energy release E g is available to GW emission, the maximum value 
of the initial amplitude is then 38 1 

1/2 



u D 



E n 



T, 



gu 



(8) 



where D is the distance of the source from the detector. Thus, the maximum strain am- 
plitude is directly determined by the energy release. However, not all the energy release is 



available for emitting gravitational waves 



7( as the energy can be dissipated by other mech- 



anisms. In fact, most of the ener gy r elease from the hadron-quark phase transition would be 



used to excite radial oscillations 39| . Moreover, the pure radial models do not generate any 
gravitational waves unless they are coupled with rotation (^J. However, because of angular 
momentum conservation, rotating neutron stars are more likely to be formed during super- 
nova explosions. If there are some other dissipating mechanisms of timescale during the 
hadron-quark phase transition in neutron stars, the fraction of energy dissipated by grav- 
itational waves is then approximately f g = 1/(1 + T guj /Td) Q. Thus, the relative damping 
timescale of the GWs with respect to that due to other mechanisms determines the fraction 
of the total energy release that can be carried away by GWs. In the literature, time scales 
of various GW modes and various energy dissipation mechanisms have been considered, see, 
e.g., refs. [3, Q] and references therein. As for many other interesting questions in the field, 
the conclusions so far are still very model dependent. 

IV. FREQUENCY AND DAMPING TIME OF THE AXIAL m-MODE OF GRAV- 
ITATIONAL WAVES 



GWs from non-radial oscillations of NSs carry important information about their internal 



structures 



41 



431 ] . One special type of normal modes of oscillations that only exists in 



general relativity is the so called w-mode associated with the perturbation of the space-time 
curvature and for which the motion of the fluid is negligible 26(. While the frequencies of 
the w-mode is significantly above the operating frequencies of existing GW detectors, it is 
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nevertheless useful to further investigate what information about the internal structure of 
NSs may be revealed from a future detection of the w-mode of GWs. 



According to Chandrasekhar & Ferrari 26| , the axial perturbation equations for a static 



neutron star can be simplified by introducing a function z(r) that is constructed from the ra- 
dial part of the perturbed axial metric components. It satisfies a Schrodinger-like differential 
equation 

g + [c 2 - V(r)]z = (9) 

where u{= u + iui) is the complex eigen-frequency of the axial w-mode and r* is the tortoise 
coordinate defined by 

where the e v and e A are the metric functions given by the line element for a static neutron 
star [26|. Inside a neutron star, the potential function V is defined by 



e 



2i/ 



V = -=-[Z(Z + l)r + 4vrr d (p(r) - P(r)) - 6m(r)] (11) 



with / the spherical harmonics index (used in describing the perturbed metric and only the 

case / = 2 is considered here), p(r) and P(r) the density and pressure, and m(r) the mass 

inside radius r, respectively. Outside the neutron star, Eq. (fTT|) reduces to 

r - 2M r / \ n 
V=— ^[Z(i + l)r-6M fl ] (12) 

where M 5 is the total gravitational mass of the neutron star. The solutions to this problem 
are subject to a set of boundary conditions (BC) constructed by Chandrasekhar & Ferrari 



26[ : regular BC at the neutron star center, continuous BC at the surface and behaving as a 



)urely 
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outgoing wave at infinity. In a recent work 25] using the continued fraction method 



451 ] . some of us have studied the frequency and damping time of the axial if- mode for NSs 
containing only the ripe/i matter. It was found that the density dependence of the symmetry 
energy has strong imprints on both the frequency and damping time of the axial w-mode. 
In this section, we examine how the appearance of hyperons and the hadron-quark phase 
transition may affect the frequency and damping time of the axial w-mode. As mentioned 
earlier, we focus on the relative effects of the symmetry energy and the bag constant on the 
first u;-mode. 

Shown in Fig. 5 is the frequency as a function of the neutron star compactness M g /R 
for the four situations considered, i.e., pure npep star, hyperonic star, and hybrid stars 

12 
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FIG. 5: (Color online) Frequency of the first to mode as a function of the neutron star compactness 
Mg/R with the symmetry energy parameter x = (left window) and x = —1 (right window). 



with two different values for the bag constant. It is interesting to see that the frequenc y is 



very sensitive to the EOS used. This is consistent with the earlier findings in refs. 25 



where a large ensemble of hadronic EOSs were used to describe the npe/i matter. Compared 
to these previous studies, it is seen that as the EOS softens when hyperons appear or the 
hadron-quark phase transition happens, the frequency increases quickly. It is also seen that 
the effect of the hadron-quark phase transition is dramatic. Comparing results in the two 
figures with x = and x = — 1, it is seen that the bag constant B plays a much stronger 
role than the symmetry energy parameter x. 
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FIG. 6: (Color online) Damping time of the first ui-mode as a function of the neutron star com- 
pactness Mg/R with the symmetry energy parameter x = (left window) and x = —1 (right 
window) . 



13 



Shown in Fig. 6 is the damping time as a function of M g /R for the four situations 
considered. It is seen that the damping time is longer for more compact NSs. Effects of the 
EOS on the damping time are appreciable but not as strong as those on the frequency. This 
might turn out to be an advantage for extracting information about the EOS of neutron 
star matter from analyzing the GW signals. 
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i i i i I i i i i I i i i i I i i i i 
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v(kHz) 

FIG. 7: (Color online) Damping time versus frequency of the first w-mode for hybrid stars. 

To assess more clearly the relative effects of the symmetry energy and bag constant, we 
present in Fig. 7 the correlation between the damping time and the frequency for hybrid 
stars. For the softer EOS with _B 1//4 = 170 MeV, the effect of the symmetry energy is small. 
On the other hand, for the stiffer EOS with B 1 ^ = 180 MeV, the density dependence of the 
symmetry energy can affect the frequency significantly. This is because the hadron-quark 
transition happens at a higher baryon density when a large bag constant is used as shown 
in Fig. [2j The symmetry energy in the hadronic phase then also plays an important role 
in determining the structure of NSs. Nevertheless, comparing frequencies with the same x 
parameter but different values of B, it is clear that the bag constant has a much stronger 
effect as it changes the underlying EOS of dense matter more significantly. 

V. SUMMARY 

Using an isospin- and momentum- dependent effective interaction for the baryon octet 
and the MIT bag model to describe, respectively, the hadronic and quark phases of neutron 
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stars, we have investigated the maximum available energy for gravitational wave emission 
due to the micro-collapse triggered by the hadron-quark phase transition in neutron stars. 
Moreover, the frequency and damping time of the first axial u>-mode of gravitational waves 
have been studied for both hadronic and hybrid NSs. Since the most uncertain part of the 
EOS of a neutron star is the density dependence of the nuclear symmetry energy in the 
neutron-rich nucleonic matter and the bag constant in the quark matter within the MIT 
bag model, we have studied effects of the symmetry energy and bag constant on the energy 
release as well as the frequency and damping time of the first axial w-mode of GWs from 
neutron stars. We have found that the energy release is much more sensitive to the bag 
constant than the density dependence of the nuclear symmetry energy. The frequency of 
the w-mode has been found to be significantly different with or without the hadron-quark 
phase transition and depends strongly on the bag constant. Effects of the symmetry energy 
and bag constant on the damping time have also been found to be appreciable but not as 
strong as those on the frequency. We have further found that the effect of the symmetry 
energy on the frequency becomes stronger with a larger value of the bag constant that 
leads to a higher hadron-quark transition density. While the predicted frequency of the w- 
mode is significantly above the bands of operating frequencies of the existing GW detectors, 
our results have indicated that the frequency of the w-mode can indeed carry important 
information about the internal structure of NSs and the properties of dense neutron-rich 
matter. 
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